arXiv:l 502.07686V 1 [math.AP] 26 Feb 2015 


THE GENERAL PEAKON-ANTIPEAKON SOLUTION EOR THE 
CAMASSA-HOLM EQUATION 


KATRIN GRUNERT AND HEDGE HOLDEN 


Abstract. We compute explicitly the peakon-antipeakon solution of the Camassa- 
Holm equation ut — utxx + 3uu.j; — 2uxUxx — uuxxx = 0 in the non-symmetric 
and Q-dissipative case. The solution experiences wave breaking in finite time, 
and the explicit solution illuminates the interplay between the various vari¬ 
ables. 


1. Introduction 
The Camassa-Holm (CH) equation 

( 1 . 1 ) Ut Utxx 3uUx ‘^UxUxx '^'^xxx — 0 

was first studied in the context of water waves in the seminal papers HIT]. It pos¬ 
sesses many interesting properties, including complete integrability and soliton-like 
solutions that interact in a manner similar to the solitons of the KdV equation. In 
the context of the CH equation, solitons go by the name of peakons and antipeakons, 
and they are the topic of interest in the current paper. The peakons, that are sta¬ 
ble solutions m, are considerably more challenging than the KdV solitons as the 
peakons experience wave breaking in finite time and become singular. 

The explicit example of the peakon-antipeakon solution for the Camassa-Holm 
equation has been a constant source of inspiration and intuition for the analysis 
of the solution of the general Cauchy problem. As in a laboratory, one can test 
one’s intuition on this particular solution that encodes most of the intricacies of the 
Cauchy problem. The key question is to analyze the behavior of the solution u near 
wave breaking where Ux{xo, t) —)■ —oo as t ^ to, yet the norm remains finite [9]. 
Multipeakons can even be used for numerical computations for the general Cauchy 
problem, see m- 

Multipeakons appear as linear combinations of single peakons of the form 

n 

u{x,t) = 

i=l 

Observe that the function u is not a smooth solution as it is not even differentiable. 
When Pi is positive, we have a peakon moving to the right, and when pi is negative 
the antipeakon moves to the left. The interesting case appears when there is at 
least one peakon and one antipeakon, which is the case of wave breaking. Explicit 
formulas exist, see, e.g., [ilZlIIlISlIlllHlIlSl- All these examples are in the so-called 
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conservative case, where the energy is preserved at the wave breaking. However, 
wave breaking allows for a dichotomy between conservative solutions and dissipative 
solutions where part of the energy is removed. The analysis of the solution near wave 
breaking requires a careful change of variables that allows for a smooth transition 
across wave breaking. For multipeakons this is discussed in [T8l[20]. Recently, a new 
class of solutions was introduced, namely so-called a-dissipative solutions that offer 
a continuous interpolation between conservative (a = 0) and dissipative (a = 1) 
solutions, see [16) . In [16] the symmetric case where the peakon and antipeakon 
completely annihilate each other at wave breaking, is analyzed in detail. 

In this paper we analyze the general case without symmetry, and where the 
solution does not vanish at wave breaking; in short we extend [161 Sect. 5] to 
the non-symmetric case. It is somewhat surprising that the non-symmetric case 
allows an explicit, albeit not simple, solution. The results are presented in this 
paper. The crux of the calculation is that one can solve exactly the equation for 
the characteristics. 

There has also been work on solitary wave solutions of the equation 


( 1 . 2 ) 


'^txx StlUx ^'^x^xx '^^XXX - O7 


which of course reduces to 0 when K = 0. The simple transformation v{t,x) = 
u{t,x — Kt) — K takes a solution u of (1.2) into a solution v of (1.1). If one wants 
solutions decaying at infinity, this transformation is of little use. However, the 
decaying solitary wave solutions of jg do not have the explicit simple form they 
have for equation ( |1.1[ ). See [HI 1221 1241 [25l [26l [27] . A complete description of 
traveling wave solutions of (1.2) can be found in [23) . 


Before we present a summary of the content of this paper, we note that the gen¬ 
eral Cauchy problem for the Camassa-Holm equation has been extensively studied 
in both the conservative and dissipative case, see niadiKniiiaiiidniiiiiii]. 
In [16] the Cauchy problem is studied in the case of a generalized Camassa-Holm 
system of the form 


'^txx “ 1 “ I^Ux 3llUx ^UxUxX '^^XXX ‘^PPx — 0 ; 

Pt + {up)x = 0. 

The present paper is organized as follows. In order to describe the a-dissipative 
peakon-antipeakon solutions, we introduce, in addition to the main unknown u, 
auxiliary variables that measure the concentration of energy. This is done in the 
form of a Radon measure p,, with the property that its absolutely continuous part 
satisfies p^c = dx. Whenever wave breaking occurs, tends to —oo and part 
of the energy is concentrated on sets of measure zero. In our case, every solution 
experiences wave breaking exactly once, and at breaking time, energy is transferred 
from the absolutely continuous part of /i to a point mass. By continuing the solution 
beyond wave breaking without manipulating p, the peakon and antipeakon are 
going to pass through each other. However, by removing an a-fraction of the 
energy that is concentrated in a point measure, the solution is either continued by a 
rescaled peakon-antipeakon for a S [0,1) or the one-peakon solution for alpha = 1. 
In addition, we introduce a measure v that keeps track of the energy changes. The 
variables (rt, p, v) are denoted the Eulerian variables. Instead of computing directly 
the solution in Eulerian coordinates, we are going to introduce Lagrangian variables 
{y,U,h,h). They are given by the characteristics y, and the Lagrangian velocity 
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defined by 

(1.3) yt=u{t,y), = u{t,y{t,^)), 

as well as two realizations of the energy given by 

(1.4) = V = y#{KOdCi- 

There are two main reasons for this change of coordinates. On the one hand this 
change of variables associates the functions h and h to the measures y and v. 
On the other hand the CH equation rewrites as a system of ordinary differential 
equations, whose solution remains smooth across wave breaking. At breaking time 
we modify this system in a continuous manner using the parameter a, as introduced 
in m- The details of the transformation between the Eulerian and Lagrangian 
variables can be found in Section and proofs can be found in [16] . In Section 
[^we compute the solution in Eulerian variables, and in Section]^ we provide the 
detailed calculations in the Lagrangian variables. However, the computation of the 
full solutions in either set of variables, requires crucial interaction between the two. 
The fact that the obtained solutions are indeed weak solutions of the CH equation 
has been established in [B Thm. 26], and is not repeated here. 

2. Transformation between Eulerian and Lagrangian variables 

The description of a-dissipative solutions of the Camassa-Holm equation is based 
on a generalized method of characteristics. More precisely, one rewrites the CH 
equation as a system of ordinary differential equations in a suitable Banach space. 
However, one faces a major problem. Energy can concentrate on sets of measure 
zero, even in the case of smooth initial data, and hence the corresponding variable 
is not a function, but rather a measure. To take care of this issue, a suitable 
change of variables needs to be introduced mapping measures to functions. Thus 
the aim of this section is to present the interplay between Eulerian and Lagrangian 
coordinates. Since these results are well-established we refer the interested reader 
to IB for the details and only state the results here. 

We start by introducing the set of Eulerian coordinates V. 

Definition 2.1 (Eulerian coordinates). The setT) is composed of all such 

that 

(i) u e iL^(K), 

(a) p is a positive finite Radon measure whose absolutely continuous part, /iac, 
satisfies 

(2.1) pac — dx, 

(Hi) V is a positive finite Radon measure such that there exists a measurable 
function f such that 

(2.2) p = fv and 0 < / < 1. 

Thus any solution of the CH equation is going to be described by a triplet 
{u{t,x), p{t,x),iy{t,x)) G V, where the measure p describes the concentration of 
energy at breaking times, while the measure v needs to be introduced for technical 
reason, but does not influence the time evolution. 

Let G be the subgroup of the group of homeomorphisms from K to K such that 

(2.3a) / — Id and f~^ — Id both belong to IT^’°°(M), 
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(2.3b) 


/{ — 1 belongs to 


where Id denotes the identity function. 

Then we can introduce the set of Lagrangian coordinates as follows. 

Definition 2.2. The set T consists of all 0 = (y, U, y^, U^, h, h) such that 


(2.4a) 

X = (C, U, Q,u^,h)e L°°{R) X [L^{R) n 

(2.4b) 

h G 

(2.4c) 

2 /j > 0, h> 0, h > 0 almost everywhere, 

(2.4d) 

lim C(C) = 0, 

(2.4e) 

y^ + h 

(2.4f) 

y^h = [/| almost everywhere, 

(2.4g) 

h > h almost everywhere. 

(2.4h) 

y + H &G 


where we denote y(£f} = f and H{t,f)= J_^h{t,f)df. 


The condition y + H G G is crucial since it in general enables to identify equiva¬ 
lence classes and hence enables to identify each element in V with one equivalence 
class in if. However, since this will not play a major role for our explicit computa¬ 
tions, we will not go into detail here. 

The change of variables between Eulerian and Lagrangian coordinates is then 
given by the following definition. 


Definition 2.3. For any inV, let 

(2.5a) y{f) = sup {y \ y)) + y <£}, 

(2.5b) h{0 = f - %(0, 

(2.5c) Ui0=uoy{0, 

(2.5d) m = foym{a 


where f is given through (2.2). Then 0 = {y,U,y^,U^,h,h) € 
by L: V —> F the mapping which to any element (u, y,,v) € T) 
{y,U,y^,U^,h,h) G F given by ([^. 


F. We denote 
associates 0 = 


For the transformation back to Eulerian variables, we apply the following defi¬ 
nition. 


Definition 2.4. Given any element 0 = {y,U,y^,U^,h,h) G F. Then we define 
(u, P-, n) as follows 

(2.6a) u{x) = U{f,) for any f such that x = y{£,), 

(2.6b) T = y#iHOdO, 

(2.6c) V = y^{h{f)dC)- 


We have that (u, y, v) belongs to V. We denote by M: F ^ V the mapping which 
to any Q in F associates the element {u,y,v) gV as given by (2.6). 
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Note that the representation in Lagrangian coordinates depends on which mea¬ 
sure V we choose. This means, given {u,^,v) and /x), then L{u,^,v) = 


{y,U,y^,U^,h,h) y^L{u,fi,y) = {y,U,y^,U^,h,h). However, if supp(/Xs) =supp(^s), 


the support of the singular measures, one can show that there exists a relabeling 


function g{^) G G such that y{0 = vigiO), U{0 = U{g{0) and h{^) = h{g{^)), 


while h{^) ^ h{g{^)). In addition, this means that the value of {u, g) only depends 
on (y, U, h), but is independent of h in this case. 


3. EULERIAN COORDINATES 


Consider the following initial data: 




(3.2a) 

(3.2b) 

(3.2c) 


qiit) = ln(L) -I- ci(t - to) - ln(ci - ‘‘’^), 

'?2(i) = - ln(T) + C2(t - to) + ln(cie^^‘“*°^ - C 2 ). 


Here to > 0 denotes the future time of wave breaking, which will take place at the 
origin. Furthermore, c\ > 0 > C 2 , and L = c\ — c^- The fully symmetric case which 


yields complete annihilation at wave breaking corresponds to ci = — C 2 [^ Note that 


qi{t) < 92(t) and 91 (to) = 92 (to)- 

Then the solution of the CH equation before wave breaking, which occurs at 
time to, is given by 

(3.3) M(t,x) =pi(t)e-l"-«i(*)l+p 2 (t)e-l"-«^(*)l, x G K, t < to. 

Define the two Radon measures by 
/x(t) = i'{t) = u^{t,x)dx 


(3.4) 



(pi(t)e -I-P 2 (t)e '^2(*^)^e^“ctx, for x < 9 i(t), 
(P 2 (t)e^“®(*^ — Pi(t)e'^i(‘)““)^(ix, for qi{t) < x < q 2 {t) 
(pi(t)e«i(*) -hp 2 (t)e«=(*)) 2 e- 2 ^dx, for q 2 {t) < x. 


Then the energy for t < to equals 
(3.5) 




9i(*) 





= 2cl + 2cl = E^. 


^Most formulas simplify considerably in the fully symmetric case, often after a limiting proce- 
dure. It is substantially easier to study the case ci = —C 2 separately. 
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Next we compute the delicate behavior at breaking time when t = Iq. The solution 
u looks like a one peakon solution with height Ci + C 2 , that is, 

(3.6) u{to,x) = (ci + C 2 )e“l“l 


as qi{to) = 92 (^ 0 ) = 0 and Pj { t ) —)■ ±00 when t —)■ 00. The special case when 
Cl + C 2 = 0 yields of course the trivial solution. 

Furthermore, we hnd that for any M C 72. measurable that 

(3.7) 

pL{t){M) = v{t){M) = f ul{t,x)dx ^ {ci + 02)^ [ - 4:CiC26o{to){M) 

Jm Jm 

where 6o{to){M) = 1 if 0 G M and zero otherwise. The part of the energy which 
has not concentrated at the origin is given by 

/ (u‘^{to,x)+ul{to,x))dx = 2 / {ci + C2Ydx 

(3.8) Jm. 7b 

= 2(ci + € 2 )^ = + 4ciC2. 


Note that 4 ciC 2 is negative and hence, comparing (3.5) and (3.8), yields that the 
amount of energy concentrated at the origin at time to equals — 4 ciC 2 . 

Thus at t = to "we find 


(3.9) 


u{to,x) = (ci + C 2 )e 

M(to-) = lim/x(to) = (ci + C 2 )^e“^l“lcta: - 4ciC2(5o(to). 
nto 


Next we introduce the parameter a G [0,1] that describes a-dissipative solutions. 
In the case of a = 0, we have the conservative solution where all energy is preserved, 
while a = 1 corresponds to the fully dissipative case where all energy concentrated 
at the origin is removed. Then an a-fraction of the energy concentrated at the 
origin at time tg is given by — 4 aciC 2 . We modify the energy concentrated at the 
origin in the measure /i while keeping it unchanged in v. More precisely, 

(3.10) v{tQ) = ul{to, x)dx - 4 ciC 2 i 5 o(to), 

(3.11) /r(to) = x)dx - 4(1 - a)ciC25o{to). 


To continue the solution for t > to we use as initial data 

(3.12) u{to,x) = (ci + C 2 )e“l“l, p{to) = (a + C 2 )^e“^l®lcix - 4 ciC 2 < 5 o(to). 


We discuss the fully dissipative case separately, and commence with the general 
case. 


The case a G [0,1): The solution will again be a two-peakon solution of the form 

(3.13) w(7,a;) =:pi(t)e-l"-«i(*)l+:P 2 (t)e"l""®(*)l, 3 . g r, t > to, 


where 

(3.14a) 

(3.14b) 

(3.14c) 


Pi(t) 


d2 - die-^^*-*°'> 
1 _ 


qi{t) = ln(L) + d 2 {t- to) 


P2{t) 


di - d2e-^^*-^o'> 
1 _ 


-ln(die-^(‘-*“)-d2), 


92(i) = - ln(T) + di(7 - to) + ln(di - ^26 ‘“^). 
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It remains to compute the values of di in terms of c^. Here L = di — c? 2 , and it 
will turn out that di > 0 > fi 2 - In particular, we have qi{t) < 92 (i) for all t > 
Furthermore, the energy for t > to equals 


(3.15) 


{t,x) + {t,x))dx = 2d\ + 2 (^ 2 . 


We are now ready to establish the connection between the pairs (ci,C 2 ) and 
(di,d 2 )- By construction we must have 

(3.16) (ci + C 2 )e“l^l = lim u{t,x) = lim u{t,x) = {di + d 2 )e~'^^\ 

t - >tQ — t ->’to + 

which implies that 

(3.17) 


dl + ^2 = Cl + C2. 


Moreover, we have for the energy, since we take out an a-fraction of the energy 
concentrated at the origin at time to, that 

(3.18) 

2dl + 2dl= lim / {u'^{t,x) + v?^{t,x))dx 

t^tQ+ 

= lim / (u^(t,a;) + w^(t,a;))dai + 4aciC2 = 2ci + 2 c 2 + 4aciC2. 

t—Ito— Jjj 

Thus, di satisfies the following quadratic equation, 

(3.19) di — di{ci + C 2 ) + (1 — q;)ciC2 = 0. 

Since we assume that di > 0, we get 


1 


(3.20) 


di = o('^i + C 2 ) + V 7(ci + C2)2 - (1 - a)ciC2 


d2 = x(ci + C 2 ) - \ -Aci + C 2 Y - (1 - a)ciC2 . 


Note that ^2 < 0, since —(1 — a)ciC 2 > 0. Furthermore, observe that in the 
conservative case with a = 0, we have dj = Cj. The new energy is given by 

( 3 . 21 ) =2{dl+dl) = +AaciC2. 

Up to this point we are only able to write down the Radon measure /i(t, x) for 
t > to, 

({pi(t)e~^^d) + for x < qi{t), 

fj.{t,x) = ul{t,x) = < {p2(t)e^~'^^d) - pAt)e^^d)-^Y(ix, for qi{t) <x < q2{t), 
y{pi{t)e^dt) + p2{t)e^^d)yg-‘^^(lx^ for q2{t) < x. 

The other measure iy(t, x) is best computed from the solution in Lagrangian coor¬ 
dinates. However, for the sake of completeness, we state the result already here 
(cf. pHit ), 

{ fJ.{t,x), foiX<qi{t), 

n{t,x) + Vm{t,x), for qi{t) < x < q2{t), 
fJ-{t,x), foiq2{t)<x. 
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where 

x ) = 4a;(l — a)cf 02(1 — 

- (di - 

It should be noted that r'm(t, [9i(i)i 92 ( 1 )]) = —4aciC2, and thus equals the amount 
of energy taken out at time t = to. Hence, the a part of the energy concentrated 
at the origin at time toj is no longer concentrated in one singular point for t > to, 
but is, in some sense, spread out over the interval [qi(t), q 2 {t)]. 


The case a = 1: In the fully dissipative case a = 1, the solution reads 

(3.23) u(t,cc) = (ci+C 2 )e-l"-("i+"^)(*-‘«^ a; e K, t > to, 

which is the one peakon solution with height and speed equal to Ci +C 2 . Moreover, 
the energy for t > tg is given by 

(3.24) / [u'^{t,x) + ul{t,x))dx = 2cl + 2cl + 4 ciC 2 = 2{ci + 02 )^■ 

Jm 


Here we can write down the measure /r(t,a:), which equals 

(3.25) ^(t,a;) = u^(t,a;) = (ci + a; S K, t > to- 


As far as the second measure 1 /( 1 , x) is concerned, we again have to postpone the 
necessary computations until the next section (see (4.29c)), but state the result 
here. 


(3.26) 


v{t,x) = n{t,x) - 4ciC25{ci+c2){t-to){t), a; e K, t > to. 


4. Lagrangian coordinates 


Next we turn to the Lagrangian variables, which for t < to are solutions of the 
following system of ordinary differential equations. 


(4.1a) 

Vt = u, 

(4.1b) 

Ut = -Q, 

(4.1c) 

II 

(4.1d) 

Ut,^ = lh + {U^-P)y^, 

(4.1e) 

ht = 2(17" - P)U^, 

(4.1f) 

hf — /li, 

where P and Q 

are given by 

(4.2) 

= \ [ e-l«(‘’«)-^(‘''')l(2t7"7/4 + /7)(t,77)dri, 

^ JR 

and 



(4.3) Q(t,^) = -1 f sign(^-? 7 )e ( 2 C/ 2 y 4 +/i)(t, 77 )^ 77 . 

4 Jr 
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The function h is an auxiliary variable whose meaning will only become clear after 
wave breaking. Note that before wave breaking, h and h coincide, that is, 

(4.4) for all ^ G K and t < to. 

However, the above system is difficult to solve directly, even in the case of the 
symmetric peakon-antipeakon solution. Instead of solving (4.1) directly, we will 
determine the solution by using the connection between Eulerian and Lagrangian 
variables directly. The key relations are 

(4.5) 2/4 = MO y, U = UOy, h = uloyy^. 

We have to determine the initial characteristic (here denoted by yo), given by 

(4.6) yo(0 = sup{y | uo{{-oo,y)) + y < ^}. 

In our example, the measure vq is absolutely continuous and hence the initial char¬ 
acteristic is implicitly given by 

rSoii) 


( 4 . 7 ) 


ujO,x)dx + yo{0 


Unfortunately, equation (4.7) is hard to solve for yoiO- However, its derivative is 

straightforward 

(4.8) 

y'oiO = ^ 


1 +^0.:r:(2/o(0) 

(1 -f (pi(0)e“'^i(°^ -|-p2(0)e“®(°^)^e^^°^^))“^, 

(1 -f (p2(0)e2^«(«)-®(o) _ p^(o)e9i(o)-y«(«))2)-i, 
(1 -f (pi(0)e«i(o) +p2(0)e'?"(°))2e-2yo(«))-i, 

(1 4 - 

(1 4- - c|e"'=2*o)2e"2SoU))-i^ 


for ^ < ^ 1 , 
for Cl < ? < 6, 
for C 2 < C 

for C < Cl, 
for Cl < C < C2, 
for C2 < C, 


where we introduced Ci and C2 as the solutions of yo(Ci) = <Zi( 0 ) and yo{^2) = 92(0), 
respectively. For reasons that will become clear later, we will benefit from having 
characteristics that satisfy yoidii^)) = 9i(0) and 1/0(92(0)) = 92(0), which is not 
automatically satisfied by ( 4 . 7 ). We use the freedom given to us by relabeling to 
modify i/o(C). To that end define 

( 4 . 9 ) f{z)=f ul{0,x)dx + z. 

J —00 

Then f{z) is a relabeling function. Observe that with this definition Ci = /(9i(0)), 
C2 = /( 92 ( 0 )) and f'{z) = u^(0, z) 4 - 1 . Introduce 


(4.10) 

which implies 

(4.11) 
and 

(4.12) 


yo{z) = yo{f{z)), 

2/0(91(0)) = 9 o(/( 9 i( 0 )) = yo(Ci) = 9 i( 0 ), 


2/0(92(0)) = 2/0(7(92(0)) = 2/o(C2) = 92(0). 
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Hence, 

(4.13) y'(e)=y'o/(0/'(e) = l. 


Thus the relabeled initial characteristic is simply yo(C) = Clearly, we could 
have chosen this function immediately, and the above argument shows that one can 
always use the identity as the initial characteristic when the initial data contains no 
singular part. However, the above argument illustrates the actual use of relabeling. 

The key observation here is that the equation for characteristics, yt = uoy, can 
be solved explicitly. Thus, the Lagrangian variables are then given, using (4.11 for 
t < to, by 


1 


for C < 91 (0), 

y(CC) = { 

1 y-mit, ^), 

for gi(0) < C < 92(0), 

\ 

[yr{t, 0 > 

for 92(0) < C, 

1 


for ^ < gi(0). 


1 y^,mit, ^), 

for (7i(0) < ^ < 92(0) 

\ 

[ y£„rit, ^), 

for 92(0) < C, 

1 


for C < 9i(0), 

tl(CC) = 

Um(t, 0 , 

for 9i(0) < ^ < 92(0) 

1 

[Ur{t, 0 , 

for 92(0) < C, 

1 

f 

for ^ < gi(0). 

h{t,^) = 1 

hm{t, 0 , 

for gi(0) < C < 92(0), 

\ 


for 92(0) < C, 


where 

yi(t,0 = C + ln(i) 

- In (l + + C2e"^‘«)e«) , 

/eC2(t-to) (ciei(‘-‘o) - C2)T>(0+ 

-ln|^ - + (cie^2(t-to) _ c2e^i{t-to))C{^) J ’ 

yr(t,C = C - ln(i) 

+ In (l + + C2e"^"‘«)e"«) , 


+ (cie - C2e + C2e°^‘‘’)e^^ L, 

y^^rnit,0 = 

X (D(0 + - C2e"i(‘-*“))C(0)”^ 

X ((cie^(*-‘«) - C2)D{0 + 
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and 


ci(i—to) —C2(t—to) 


L + (cie *o) _ cj^gcito _ 226 *o) + c2e'^2*o)e? ’ 

Um{t,i) = -c2) + 2C(0^?(0^"e"^^‘-*“Hci +C2) 

X ((cie^(*-‘«) - 02)15(0 + L^e"i(‘-*'>)C'(0)"^ 

X (i5(0 + - C2e"i(‘-*“))C(0)”\ 


C/,(i,0 = 


^2gCi(t-fo) _ g2gC2(t-fo) 


and 


L + (cie^l(* *o) _ Cj^g Clto _ g 2 gC 2 (t to) _|_ g^g C2fojg 4 

(L + + c2e'^2*o)g5)3 ’ 

^m(^,0 = U‘^it,Oyd't^O - 4pi(0P2(i)e‘^^^‘^"'^"^‘0«(^>C)> 

(g 2 gCi(t-fo) _ g2gC2(t-to))22^g-25 

4 j ~ 




hi{t,C) = 


(L + (Ole'll*o) — Cie ^1*0 _ g2gC2(t to) _|_ g2g coto^g 5^3 


Here 

(4.14) 


C(0=C2-cie-^‘“+Lef+"^‘«, 

£>(0 = - Lcie« + Lc2e^-^*°. 

The above formulas are obtained as follows. One starts by computing the charac¬ 
teristics using (3.3). The cases ^ < gi(0) and 92 ( 0 ) < C are more or less straightfor¬ 
ward. However, the challenge is to solve yt{t,C) = u{t,y{t,^)) for gi(0) < ^ < (^ 2 ( 0 ), 
since one has to rewrite the resulting ordinary differential equation in a suitable 
way by applying several changes of variables. The equation reads 

L 




r(e 


1 _ gL(t-to) 

Introducing the equation becomes 


ci(f-to)-y(t,5) _ g-c 2 (t-fo)+y(f,?) 


)• 


xtit,0 = 




Cl(t-to) _ g-C 2 (t-to ).^2 




1 _ gL(t-to) 

By further introducing vit,^) = -L{1 - we find 

= -R'it) + Rit)vitA) + 

where 

R(A - - C2 

^ > 1 _ gi(i-to) ■ 

Defining the logarithmic derivative v(t,^) = —St(t, ^)/s(t, ^), the equation reads 
sttit,0 = {R(.'t)s{t,C))t, or stit,^) = Rit)s{t,0 + I{^). 

Yet another substitution = (1 — turns the equation into 

= -C2w{t,^) + - e^(*-‘o)) 
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where /(^) is a constant of integration. Finally, by defining z{t, 
we obtain 

which can easily be integrated. By returning to the original variables, we find 
what is denoted If is not clear at first sight that the characteristics for 

91 ( 0 ) < C < 92 ( 0 ) are well-defined, and one has to check that the argument in 
the logarithm is positive and bounded. In particular, one can show that D{^) + 
(cieC 2 (t-to) _ < 0 and - C 2 )D{^) + < 0 

for all ^ G [gi(0),(72(0)]. 

Once the characteristics are known for t < to, the Lagrangian velocity 

= u{t,y{t,^)) is straightforward. 

As far as the energy variable h{t,^) is concerned, the necessary computations 
simplify considerably by rewriting the equations. For ^ < <71 (0), one has that 
Ux{t,yit,^j) = -u{t,y{t,^)), and hence h{t,^) = Similarly, for 

92(0) < d one has that ■Ua;(t, y(t, ^)) = u{t,y{t,^j), and hence /i(t, 5) = U^{t,^)y^{t,^). 
Again the challenging case is (7i(0) < ^ < 92 ( 0 ). First we calculate the limit as 
t —)■ to—- Thus it sufhces for our purposes to observe u‘^{t,y{t,^)) = v?{t,y{t,d)) — 

and hence/i(t, = C/^(t, 0- 

The representation we have chosen so far makes it quite easy to compute the 
limits as t ^ to- 


( 4 . 15 a) 

yihd) 

( 4 . 15 b) 

Uito,0 

( 4 . 15 c) 

h{to,0 


C + ln(T) 

— ln(L -I- (ci — — C2 + C2e’^^*°)e^) ’ 

< 0 , 

e - ln(T) 

+ ln(L -I- (ci — — C2 + C2e~‘^^*°)e~^) ’ 


for ^ < (7i(0), 

for gi( 0 ) < ^ < (72(0), 

for (72(0) < 


(ci + C2)e^(*'>’^), 
Cl -b C2, 

(ci + C2)e~y^*°'^\ 


for ^ < gi(0), 

for gi(0) < C < <72(0), 

for 92(0) < d 


' _ (ci+C2)^Z/^e^^ _ 

(L+(ci —cie*^i *^0 —C2+C2e‘^2f^0)e^’ 

_ 4c^c^(l-e~^*o)^e^ _ 

(Le“'^i*o +C2 —cie“^*o+(—ci+C2e“^*o+Le'^2fo )e^)2 ’ 

_ {ci+C2fL’^e-^( _ 

(L+(ci —cie“^i^o —c2+C2e“^2fo)e“ 5)3 5 


for ^ < gi( 0 ), 

for gi( 0 ) < ^ < 92(0), 

for 92(0) < d 


All these limits are obtained by direct computations. As far as h{to ,0 fo'' ? € 
[(71 ( 0 ), (72(0)] is concerned, observe that y^{t,^) —> 0 for all ^ G [(7i(0), (72(0)] as t ^ 
to-- Thus, limt^tn h{t,^) = - 41 imt^toPi(t)p 2 (t)e«i(*)“®(*L?(f>^)> since U{t,^) G 
i?^(IR) for all times, which simplihes the calculations considerably. Moreover, note 

that h{to ,0 = (ci +C2)^e"2|ydo,C)ly^(7g^^) for ^ ^ [gi(0),92(0)]. 

Next we are going to show that also in Lagrangian coordinates we have that the 
amount of energy concentrated at the origin at time t = to equals —4ciC2. Note that 
G M I y{to ,0 = 0 } = [<71 ( 0 ), (72(0)], and thus the amount of energy concentrated 
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at the origin is given by the following integral 


ri2{o) 


' 9 l( 0 ) 


h{to,Od^ = 




— Cl + C2e~^^° + Le= 2 ‘o / cic2(l-e--^*0)2 


ciC2e-‘^2tO(l- 


to 


jjZ 


= -4ciC2. 


Let us now compute the solution for t > to: 


The case a = 1: The case a = 1 yields the dissipative solution, thus we will 
introduce the function h{t, for t > to as follows 


(4.16) 


h{t,^) = 


for gi(0) < f < 52(0), 


h{t,^), otherwise. 


In addition, we have to redefine our system (4.1) of ordinary differential equations. 


To be more explicit, we have to replace h{t,^) by h{t^^) everywhere on the right- 
hand side for all t > to- Note that this also means that h{t,^) has to be replaced 
by h{t,^) in ( |4.2[ ) and ( |4.3| ). Observe that even if h{t,^) differs from h{t,^) only 
on the interval [51 (0), 52(0)], the solution will be influenced for all ^ G K, due to 
P{t,^) and Q{t,^) on the right-hand side of (4.1). We now want to conclude that 
the solution for f > tg is given by a one-peakon traveling wave with height ci -I- C2. 
As far as our system of ordinary differential equations is concerned, it reduces, for 
C G (9i(0)>92(0)) and t > to, to 


(4.17a) 

yt = u. 

(4.17b) 

Ut = -Q, 

(4.17c) 

yt,i = 0, 

(4.17d) 

= 0, 

(4.17e) 

ht = 0, 

(4.17f) 

ht = 0. 


Thus y({t,^) = 0 , U^{t,^) = 0 , and h{t,^) =0 for all ^ G [51(0), 52(6)] and t> to- In 
particular, y{t,^) = c{t), U{t,^) = c'(t), and Q(t,^) = c"{t) for all ^ G [51(0), 52(0)] 
and t > to, where c{t) denotes some suitable function only depending on t. Since 
both y{t, •) and U{t, •) are continuous for any t > to, we have 


(4.18) y{t,qi{0)-) = y{t,q2{0)+) and U{t,qi{0)-) = U{t,q2{0)+), t>to. 


As far as the solution for ^ ^ [51(0), 52(6)] is concerned, we are using a squeezing 
and relabeling argument based on the considerations for ^ G [51(0), 52(0)] so far. 
Namely, let for t > to, 


(4.19a) 

(4.19b) 

(4.19c) 


y{t,0 

uit,0 

fi{t,0 


y{t,0, forC<5i(0), 

y{t,^ +92(0)-qi{ 0 )), for 5 i( 0 )<^. 


U{t,0, forC<9i(0), 

U{t,^ + q 2 { 0 )-qi{ 0 )), for 5 i( 0 )<^, 


h{t,^), for^< 5 i( 0 ), 

^(L'C+ 92 ( 0 ) - 5i(0)), for5i(0)<^, 
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(4.19d) 

(4.19e) 

(4.19f) 




Pit,0 = 


Q{t,0 = 


P{t,^ + q2{0)-qim, 

Q{t,^ + q2{0) - <7i(o)), 


for ^ < gi(0), 
for 9 i( 0 ) < C, 

for ^ < 9 i( 0 ), 
for gi(0) < 


which means that we have taken out the part of the function where the energy is 


concentrated. However, due to (4.17) and (4.18), both and 


Q{t,^) are continuous. In particular, the triple {y{t, ^), h{t, ^)) satisfies the 

following system of ordinary differential equations for t > Iq, 


(4.20a) 

(4.20b) 

(4.20c) 

(4.20d) 

(4.20e) 

(4.20f) 


yt = U, 

Ut = -Q, 

yt,i = 

Ut,i = \h + {U^-P)y^, 
ht = 2{U^ - P)U^, 
ht = ht- 


A close look reveals that the above system coincides with the one describing the 
conservative solutions of the Camassa-Holm equation. To be sure that the solution 


of (4.20) for t > to coincides with the one-peakon solution with height ci +C 2 , there 


are two more properties we have to check. On the one hand, it is left to show that 
(4.21) f{0 = y{to,0+ [ Hto,v)dr] 

is a relabeling function, which guarantees that belongs 

to P, the set of Lagrangian coordinates. On the other hand, if we can also check 
that y{to,^) is a relabeling function, we can map y{tQ,^) to the identity, thereby 
ensuring that (y(toj C)i C); ^(^o, 0) is a Lagrangian representation of the one- 
peakon centered at the origin. Both claims follow from applying [M Lemma 3.5]. 
This means, in particular, according to ( |3.23| ), that 


C 2 )' 


,-\v(t,i)-(ci+C2){t-to)\ 


^ e K, t>to, 


(4.22) U{t, p) = u{t, y{t, 0) = (ci 
and especially 

= sgn{y{t,^) - (ci +C 2 ){t - to))U{t,0{U{t,^) -a- C 2 ). 

Let z{t,^) = y{t,^) - {ci + C 2 ){t - to) and V{t,^) = U{t,^) - (ci -I-C 2 ). Then the 
first two equations in (4.20) rewrite as (see Figure 

(4.23a) Zi = V, 

(4.23b) Vt = -sgn(z)V(V + (ci+C 2 )), 

and the above system has a unique solution in W ^’°°(]R) X IF^’°°(M). This means 
in particular that we can solve the above system for every ^ £ K. 
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Figure 1. The vector field for the functions (z(t),V(t)) in (4.23) 
for Cl + C 2 = 1 . 


Given ^ S K such that ^ (0,0), then ^ (0,0) 

for any finite time t > tg. In particular, if {z{tQ,^),V{tQ,^)) = (0,0) for some 
^ e K, then {z{t,^),V{t,^)) = (0,0) for all t > to- This means, in particular, that 
the peakon is traveling along a characteristic, i.e.. 


(4.24) y(t,gi(0)) = (ci + C2)(t-to), 

(4.25) G(t,(?i(0)) = (ci+C 2 ). 


Furthermore, we have for all ^ S [9i(0), 52 ( 0 )], due to (4.18|) and (4.19), that 


(4.26) y{t,0 = 2/(^>9i(0)-) = y{t,qi{0)) = (ci + C 2 ){t - to), 

(4.27) Uit,0 = C/(t,( 7 i( 0 )-) = Uit,q,i0)) = (ci +C 2 ), 

(4.28) Qit, 0 = Q(t, qi (0)-) = Qit, qi (0)) = 0. 


As a byproduct of our analysis we obtained a lot more information about our 
solution in Eulerian coordinates. Indeed, we have for t > to, 


(4.29a) uit, x) = (ci + , 

(4.29b) fi{t,x) = ul,{t,x)dx, 

(4.29c) iz{t,x) = M^(t,ai)da; - 4ciC2(5(ci+c2)(t-to)(i)> 


since the peak is traveling at speed Ci + C 2 . 

The solution in Lagrangian coordinates for t > to reads. 


(4.30) 

? + ln(T) 

- ln(L + C 2 e="*«)e«), 


y{t,0 = (ci + C2)(t-to), 


e-ln(L) 

+ ln(L + (Le('=i+^ 2 )(i-to) _ cie-^i‘“ + C 2 e-="*“)e-^), 


for ^ < gi(0), 

for qi(0) < C < 92 ( 0 ), 

for 92 ( 0 ) < C, 
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(4.31) 

uit,0 

(4.32) 

(4.33) 


(ci+C2)Le-(=i+=2)(‘-*o)e« 

L+(Le-(‘'i+“2)(‘-‘o)-cie‘=i‘o+c2e‘=2‘o)e5 ’ 

< (Cl + C 2 ), 

(ci+C2)Le(‘=i+‘=2)(‘-*o)e-« 

^ L+(LeC'=i+‘=2)(t-‘o)-cie“‘'i‘o+c2e-‘=2‘o)e-« ’ 


for ^ < qi{0), 

for gi(0) < ^ < (? 2 ( 0 ), 

for 92 ( 0 ) < C, 


' (ci+c2)^L^e-^<‘^i+“2)(*-*o)e2« 

(L+(Le-(“i+‘'2)(‘-*o)-cie‘=i‘o+c2e‘=2‘o)e«)3 ’ 

4c^ c|(l — 

(Le“^i*o+c2 —cie“^*o +( —ci+C2e“^*o +Le‘^2*o )e€)2 ’ 
_ (ci+C2)^L^e^<‘^i+°2)(*-*o)e~^^ _ 

. (L+(ie(‘=i+‘=2)(*-‘o)-cie-‘=i‘o+c2e“‘=2‘o)e-£)3 ’ 


for ^ < gi(0), 

for gi(0) < ^ < 52 ( 0 ), 

for 92 ( 0 ) < 


_(ei+C2)^£^e~^<‘^i+°2)(*-*o)e^£_ 

(L+(Le-(=i+‘'2)(‘-‘o)-cie‘=i‘o+c2e‘=2‘o)e£)3 ’ 

0 , 

(ci+C2)^L=*e^t‘^i+°2)(*-*o)e-^£ 


for ^ < gi(0), 

for gi(0) < ^ < 92 ( 0 ), 


I.C1+C2J i. e ' -1 ' u.e ’ fnr « f'n'l <' .C 

^ (L+(ie(‘=i+‘=2)(*-‘o)-cie-'=i‘o+c2e-‘=2‘o)e-£)3 > ^2 V '27 ^ • 


Next we study the general case. 


The case a G [0,1): Here we have to introduce the function h{t,^) for t > to, as 
follows 


(4.34) 


h{t,^) = 


h{t,^) - ah{to,0, for gi( 0 ) < ^ < 172 ( 0 ), 
otherwise. 


Observe that in the fully conservative case, a = 0, we have h{t,^) = We 


have to redefine our system (4.1) of ordinary differential equations for t > to. To 


be more explicit, we have to replace h{t,^) by h{t,^) everywhere on the right-hand 
side for all t > tg. Note that this also means that the function h{t,^) has to be 
replaced by h{t, ^) in the definitions of P and Q, i.e., in ( |4.2[ ) and ( |4.3| ), respectively. 
Note that for a G (0,1), h{t, ^) differs from h{t, only on the interval [( 71 ( 0 ), ( 72 ( 0 )]. 
However, due to the nonlocal nature of P{t,^) and Q{t,^) on the right-hand side 
of (4.1 1 , the solution will be influenced for all S M! 


Again it is difficult to solve (4.1) with h{t, replaced by h{t,^) on the right-hand 


side for t > to- Thus we proceed as follows. We know that the solution in Eulerian 
coordinates at breaking time is given by 

(4.35a) u{to, x) = {di + ^ 2 ) 6 “'“^' = (ci + C 2 )e“l'^l, 

(4.35b) 7 £(to, a;) = M^(to, x)dx - 4 did 2 < 5 o(^o) = x)dx - 4(1 - Q;)ciC 2 do(io), 
(4.35c) v{to,x) = u^{to, x)dx — 4 ciC 2 do(to)- 


In the special case a = 0, we obtain the well-studied conservative solution, that 
j for j = 1, 2 and the solution for t > to equals (3.13) with pj = pj 


IS, dj = c 

and qj = qj for j = 1, 2. In other words, the solution (3.3) is valid for all t > to- 
Moreover, in this case h{t,^) = h{t,^) for all t > to and hence p{t) = v{t) for all 
t > to- Thus the measure v is not needed as it does not add any information. 
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For a G (0,1) the situation is a bit more involved. Let U{t, ^), h{t, ^), h{t, ^)) 

be the a-dissipative solution in Lagrangian coordinates. Since we replaced h(t, 
by h{t,^) everywhere on the right-hand side of (4.1) for t > to, we have that the 
time evolution of is independent of Thus if 


(4.36) 


yito,0 


h(to, y)dr] = y(to,^) + H{to, e G, 


then {y(t,^), U(t, ^), h{t, ^), h{t, ^)) is the solution of ( |4.1| ) with (valid!) initial data 
iyito,0,U{to,O^Hto,O^Hto/oUo^t> to- However, {y(t,^),U(t,^),h{t,^),h{t,^)) 
in Lagrangian coordinates corresponds to the conservative solution in Eulerian co¬ 
ordinates with initial data {u{to), yL{to)), which is given according to the case 
a = 0, by (3.3) with Cj replaced by dj for j = 1, 2. 

Thus it is left to show that y{tQ, -G H{to, ^) is a relabeling function. We apply 
the fundamental lemma [161 Lemma 3.5], which reduces this difficult task to showing 
that there exists c > 0 (which may depend on fg) such that 

(4.37) c < y^{to,^)+ h{to,^), for all ^ S M. 

To that end we observe that 

W„.„ciio „^„c 2 to\ for^<gi( 0 ), 




(4.38) 


and 


ydto,0 > { 0 , 




for 9 i( 0 ) < C < 92 ( 0 ), 
for 92 ( 0 ) < C, 


(4.39) h{to,0> 


Thus 


0 , 

4did2e'^^^^^ min(ci — C 2 e“^*°, 

0 , 


i;(ci 


pCito _ 


C2e 


C2to 


), 


ydtod)+Hto,0 > 4 did 2 e«i(°) min(ci - C 2 e 

j;(cie"'=i‘o - C2e"'=2*“), 


for ^ < qi( 0 ), 

for gi( 0 ) < C < 92 ( 0 ), 

for 52 ( 0 ) < C- 

for ^ < gi( 0 ), 
for gi( 0 ) < C < 92 ( 0 ), 
for 92 ( 0 ) < 


which proves (4.37), since all terms on the right-hand side are strictly positive due 
to our assumption that Ci > 0 > C 2 and hence also di > 0 > ^ 2 - 

Finally recall that = h{t,^)+ah{to, ^), which at first enables us to compute 

h{t,^) and in a further step to derive v{t). 

We find that the solution in Lagrangian coordinates reads 


1 

^yidd), 

for ^ < 9i(0), 

ydd) = { 

1 2/m(L C)j 

for 9i(0) < C < 92(0), 

\ 

[yr{t, 0 , 

for 92(0) < C, 

1 


for ^ < 9i(0), 

ydtd) = { 

1 y^,mit, ^), 

for 9i(0) < ^ < 92(0) 

\ 

[yi,rit, 0 , 

for 92(0) < 

1 

f 

for ^ < 9i(0), 


Umdd), 

for 9i(0) < C < 92(0), 

1 

[Udtd): 

for 92(0) < 
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Here 

yr{t,0 

and 

yuit^O 

yi,r{t,0 

and 

Ui{t,0 

Um{t,0 


Ur{t,0 

and 



for ^ < 9i(0), 

L(t,0 = < hmit,0, 

for gi(0) < ^ < 92(0), 

[hr-{t,^), 

for 52(0) < 

ihi{t,0, 

for ^ < 9i(0), 

h{t,Cl = < hmit,0, 

for gi(0) < ^ < 92(0), 

[hr{t,0, 

for 52(0) < 


5 + In(LZ) 

- \n{LL + (Z(-cie"i‘« + C2e==^*“) + 

In ( - - 1) + (-^i + d,e-^(^-*o^){S{0 + 1) \ 

^ - \n{LL) 

+ \n{LL + (L(-cie-‘=i‘« + 026-==^*“) + - d2e''^(*-‘«)))e-«), 


LL 

LL + (L(—Ole'll *0 + C2e'^2‘o) + — (i2e“‘^^^‘“*“^))e?)' 

-2did2e^^^*-^°\l - 

X (l(5(0 + 1) + - die‘^^^^-*o'>){S{0 - 1))”' 

X ((di - d2e-^^^-*«^){S{0 + 1) - Le‘^^(*-‘»)(5(e) - 1))”', 

LL 

LL + (L(—cie“'^i*o + C2e“'^2to) -)_ L{die‘^^^*~*°^ — d2e‘^'^^*~*°'>))e~^ ' 

-d2e-rf2(i-io))i^eC 

LL + (L(—Ole'll*0 + C2e'^2to) _(- L(die“‘^i(‘“*o) — (i2e“‘^^^‘“*“^))e? ' 

L((d2 _ d2e-i(t-io))(5(^) + 1)2 _ 2(^2 _ rf2)gd2(t-to)(^(^)2 _ i) 

+ e‘^^(‘-*»)(d2e''2(t-to) _ d2gdi(t-to))(5^(^) _ 1)2^ 

X ((L(S'(C) + 1) + - die'^=^(*-*'>))(5'(C) - 1))”^ 

X ((di - d2e-^(*-‘«))(5(e) + 1) - Le‘'^(*-‘»)(5(e) - !)))”', 

(d2gdi(i-to) _d2e^2(*-‘«))Le-« 

LL + (L(—cie“'^i*o + C2e“'^2to) -|- L{dieLd-to) — d2e‘^2(t-to)))e“5 ’ 

- die-'^"(*-‘«))2L3Le2« 

[LL + (L(—cie^i*° + C2e'^2to) _(- L(die“‘^i(‘“*o) — d2e“'^2(t-to)))e4)3 ’ 
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hr{t:0 


UHt,0ydt,0 - 

(^2g<ii(t-to) _ ^2gd2(i-to))2^3j^g-2C 

[LL + (L(—Cie“°i*o + C2e~‘^^*°) + ’ 


and 


hi{t,0 

hm{t,0 

hr{t,0 


[LL + (L(— 016 *^ 1*0 + C2e'^2‘o) + ’ 

U‘^(t,0yd't^0 - - 2aciC25''(C), 

(^2gdi(t-to) _ ^2gd2(t-to))2^3|^g-2C 

(LL + (L(—cie“^i*° + C2e“'^2to) -|- L((iie‘^i(*“*t>) — (i2e‘^2(t-to)))e“5)^ ’ 


where 


Q(^o+,C) = (1 - a)Q(io-,0 

'2ciC 2L(1 - e-^*'')2 - (ci - C2e-^‘« + Le‘=2‘«)(LC'(0 + £>(C)) 


■ 


(—Ci+C 2 e + Le‘^^*o){LC{^) + D{^)) 


and 


sd) = 


Q{to+, C) 2ciC2L(l - - (ci - C2e-^*« + Le'=^*o)(LC(C) + D{ 0 ) 

did2 (-Cl + C2e-^*o + Le^2to)(LC(0 + ^(0) 

Note that S'(^) e [—1,1] for ^ S [gi(0), 92 ( 0 )]. Moreover, direct computations yield 

L2(1 - 


(4.40) 


5'(e) = -2 ciC2 


{LCiO + DiO?- 


(4.41) 


In this case the solution in Eulerian coordinates reads for t > to (recall (3.13)) 
u(t,x) =:Pi(t)e-l"-«i(‘)l +P2(0e”'"”®^‘^', 

n{t,x) = u^{t,x)dx, 

ul{t,x)dx, ior X < qi{t), 

v{t, cc) = .^ M^(t, x)dx + Vm{t, x)dx, for qi{t) < x < g 2 (i), 
u‘l{t,x)dx, for q 2 {t) < X, 


where 

Vm{tjx) = 4 a(l — a)ciC2(l — 

^ j'gX-di(t-to)(|^ _ rf^gd2(*-to)) 

- (di - d2e-^(‘-*“) - Le'^^^*-*")) 

The solution u is displayed in Figure!^ while U is plotted in Figure]^ The char¬ 
acteristics y are visualized in Figure and the measures can be found in Figure 

El 
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Figure 2. The function u is plotted for times t = —1.5, t = 1.0 = 
to, and t = 3.0. Parameter values ci = 0.8, C 2 = —2.0, a = 0.5. 
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